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The multiphoton resonance transitions between ground hyperfine states are used for the time 
modulation of the electron capture decay of hydrogen like ions with the Gamow- Teller transition 
1 + — 5- + . The proposed mechanism offers a time oscillating decay with the frequency up to 0.1 
Hz. The experiment to observe the modulation is proposed for ions stored in a Penning trap. An 
attempt to understand the GSI anomaly with multiple photon transitions is made. 
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Longstanding attempts have been made to change ar- 
tificially the decay of radioactive nuclei using different 
external factors. It was suggested by Segre [l| that the 
orbital electron capture (EC) decay rates depend on the 
density of atomic electrons within a nucleus. Thus, en- 
vironmental conditions may alter electron densities and 
affect EC decay rates. However this effect is small and 
does not exceed the fractions of percentage point. 

The introduction of Penning traps and storage rings 
into nuclear physics has allowed to study radioactive 
hydrogen-like (H-like) ions. Recently the single ion spec- 
troscopy technique with time resolution less than 1 s has 
been developed at GSI Darmstadt. With this elegant 
method, the oscillatory modulation of the exponential 
EC decay for the 14 0p r 58+ and 142 Pm 60+ ions was found 
in 0. The modulation period of 7.1 s is incommensurable 
with all the energy intervals involved in the experiment 
with exclusion of the neutrino mass difference. Perhaps 
this is the main motivation of the Letter [3| which relates 
the anomaly to neutrino mixing. 

Our motivation for this work is to show that there 
is a mechanism related to Rabi oscillations which offers 
such a small period. To observe the time modulation of 
the EC decay, we propose to apply Penning trap mass 
spectrometry with ions in the intermediate mass region 
having the "sterile" relative to the EC decay hyperfine 
state. We analyze also the effect of the magnetic fields 
of dipole and quadrupole magnets in the Experimental 
Storage Ring (ESR) on the GSI oscillations. 

The lsi/2 state of a H-like ion has two hyperfine levels 
with the total angular momentum F — I ± 1/2, where 
I is the nuclear spin. In the case of a positive nuclear 
magnetic moment /z, the ground state of an ion has F = 
7 — 1/2, whereas for negative /i, the order of these energy 
levels is reversed. When 1 = 1, the hyperfine splitting of 
the lsi/2 state is given by [1] 
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where m, M, /^jv, and a are the masses of electron and 
proton, nuclear magneton, and the fine-structure con- 
stant, respectively. For nuclei with I* = 1 + , the EC 
decay from the F = 3/2 state is forbidden because fully 
ionized daughter nuclei have I* = + . This feature of 
such ions may be used for observing the time modula- 
tion of the EC decay. 

The oscillatory behavior of the EC decay of such ions 
can be obtained with the Rabi resonance method used in 
nuclear magnetic spectroscopy. We consider an ion in a 
static magnetic field in the z direction, resulting in the 
Zeeman splitting of both hyperfine levels. Transitions 
between the Zeeman components of these two states are 
driven by the transverse oscillating field along the x direc- 
tion with frequency uj. When the frequency u is tuned 
to AE/%, the oscillation of populations for the states 
F = 1/2, M and F = 3/2, — M, where M = 1/2 or -1/2, 
is modulated with the Rabi frequency, which is propor- 
tional to the amplitude of the driving field. Because the 
EC decay probability is proportional to the occupation 
of the F = 1/2 state, the Rabi frequency may be ob- 
served as the time modulation of the EC decay rates. 
The time modulation of EC-decay is also possible with 
the resonance transitions involving several field quanta. 

Penning trap. — A single H-like ion in a Penning trap 
is an ideal object for observing the modulation of the 
EC decay. The ion having a charge-to-mass ratio q/M is 
confined in a strong magnetic field B(0, 0,Bq) superim- 
posed with a weak electrostatic quadrupole field [5[ . The 
solution of the motion equations yields the cyclotron, ax- 
ial and magnetron harmonic oscillations. The ion mass 
is determined via the determination of the free-space cy- 
clotron frequency uj c = qB z /Mc which derived from three 
measured harmonic frequencies. The single ion spec- 
troscopy developed in the ESR experiments can be in- 
stalled at a Penning trap fgj . In this technique, it is pos- 
sible to realize the fully controllable experiment on the 
time modulation of the EC decay by single- and multi- 
photon resonance. 

We are interested in the time evolution of the two hy- 
perfine states of an ion exposed by the trap magnetic 
field B and irradiated by the in-plane elliptically polar- 
ized driving field h(b x sin uit, 0, b z cos uit). The strength 
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of the latter is restricted by the frequency resolution of 
the single ion spectroscopy. It is necessary for a shift in 
the cyclotron frequency w c caused by the influence of the 
field b to be small compared to the change in uj c due to 
the transition from a parent ion to a daughter nucleus. 
As a rough approximation, this gives b < BqQec/(Mc 2 ), 
where Qec is the decay energy. 

The total Hamiltonian is H = H + U(t) + V[t), 
where H = Hp + (e/2)[ar]-B describes the lsi/ 2 state 
of a H-like ion in the Penning trap. Here Hp is the 
Dirac Hamiltonian with hypcrfine interaction and a is 
the Dirac matrix. The eigenfunctions and eigenvalues of 
H are |+) = \3/2,M + ), E 3/2 + {l/3)gn B B M+ = Huj 3 
and |-) = \1/2,M), E 1/2 - (l/3)g(j, B B M = hu>i. Here 
fiB is the Bohr magneton and g is the gyromagnetic ra- 
tio of the bound electron in a H-like ion. The terms 
U(t) = (e/2)[otr] z b z cos uit and V(t) = (e/2)[ar] x b x sinwi 
represent the diagonal and off-diagonal parts of the in- 
teraction between an ion and the driving field. 

There are three independent transitions between the 
lower and upper Zeeman sub-states which can be used 
for the modulation of the EC decay: 

(i) The transitions between the states 1 3/2, +1/2) and 
|l/2, —1/2) or |3/2,-l/2) and |l/2,+l/2) with the en- 
ergy AE = E 3 / 2 — Ei/2 = huio regardless of the sign of fx. 

(ii) The transitions between the states |3/2,+3/2) and 
|l/2, +1/2) for (i > (upper sign) and jU < (lower 
sign) with the energy AE — (2/3)g^B-Bo- 

(iii) The transitions between the states 1 3/2, ±3/2) and 
1 1/2, ±1/2) for \i > (upper sign) and [i < (lower sign) 
with the energy AE + (2/3)<7/Zb-Bo- 

Thus, we have the two discrete states |+) and |— ) with 
the energies e 3 = hui 3 and ej = h(uji — iX/2), where A is 
the EC decay constant. Let the amplitudes for an ion to 
be in these states are a+(t) and a_(i). Then the system 
evolves according to the equations 
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To eliminate the diagonal oscillating terms in ([2]), we 
perform the transformation 

a± = A± cxp [— i(u>i + uj 3 )t/2 — i(U±/fuj) smut], (4) 

which reduces the equations for the transitions (i) to the 
form 

i A- = - \ (lo + i A) A_ + VA + sin ut, 
iA + = \uj A + + VA- sinwi, (5) 
The two other transitions are described by the equations 
iA- = — ^(CjQ+i\)A-+VA + exp (i^ sinwi) sinwt, 



iA + = 7}UJqA + + VA- exp smut) sinwt, 
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where ujq = u> 3 — u)\, U = (U- — U+)/h. Through the use 
of the expansion of an exponent in the Bessel functions, 
these equations are converted to 
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where J' m is the derivative of the Bessel function. The 
equations (JS|) and ([7]) describe an ion with a positive nu- 
clear magnetic moment fi. When ji < 0, it is necessary 
to change the sign of loq and luq. Up to this point all of 
the equations are exact. 

The equations (0 for A — are the basic ones in the 
theory of the resonance interaction of the two level sys- 
tem with a periodic field. When u ~ ujq, we get the usual 
Rabi formula for amplitudes replacing the sine by one ex- 
ponent [3- For the K-photon resonance when cjo ~ Kui, 
Shirley [8| found the approximate solution for a weak 
driving field. In the adiabatic approximation, K 1, 
the solution of §5§ was first obtained by Zaretskii and 
Krainov [9| and then by Duvall et al. [lQ]. In the weak 
field regime, b x <ti Bq , which is adequate for the Penning 
trap experiment, we use Sambe's perturbation approach 
[111 ] provided that In the K-th order of the per- 

turbation theory in the quasienergy representation, Eqs. 
([5]) can be approximated by 
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is the shift of the resonance frequency related to the 
dynamic Zeeman effect which is referred as the Bloch- 
Siegert shift. The shift is of the order of coo(b x /Bo). For 
K = 1 the coefficient before round brackets is equal to 
1/4. The half width of if-photon resonance is 
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the photon number K being odd for these transitions. 

We now turn to a close examination of Eqs. (J7J for the 
transitions (ii) and (iii). When the frequency Kuj (K is 
even or odd integer) is nearly resonant with the energy 
separation luq of the two states, the term with m = K 
represent a resonance driving field, whereas those with 
m / K result in the resonance frequency being shifted. 
In the second order perturbation theory, we obtain 
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Substitution of A± by A± exp[+i(Ku)t — 7r/2)/2] leads us 
to © with co res given by (jlip and the half width 



T K 



2g[i B b x 



J'k 



( 2g\i B b z 
V 3hco 



(12) 



We assume that an ion is in the lower hyperfine state 
when it is injected in a Penning trap at t = 0. With this 
initial condition Eqs. (|SJ) can be easily solved for A <C Tk- 
The probability of finding an ion in the decaying state |— ) 
is given by 
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where 57 = y/ A 2 + r 2 ^ is the Rabi frequency and A = 
w. res — -Kw is the detuning. The probability depends on 
the lower state angular momentum F — F± . Expression 
(fT3| describes damped oscillations with the frequency f2 
decaying on the scale of half the EC decay constant A. 
Nevertheless, the total number of decays for infinite time 
calculated with (fT3|) is equal to the number of radioactive 
ions at t = 0. 

To elucidate the time modulation of EC decay, we uti- 
lize the concept of the dressed atom introduced by Shirley 
Q as the quantum limit of the semiclassical theory. The 
ion with \i > coasted in a trap and exposed to a reso- 
nance field is described by the quantum wave function 

*(t) = a+(t)\+)\N-K) + a-(t)\-)\N), (14) 

where \N) and \N—K) are the wave functions with the 
fixed number of photons having the energy Hlu. Thus, 
the dressed ion has the two hyperfine states entangled 
with photons. The EC decay at the time t disentangles 
the state (fL4")l by projecting ^(t) on the state |— ) with 
the probability Pp_ (t) ([T3"| , to which the decay rate is 
proportional. 

To observe the time modulation, the H-likc radioactive 
ion has to have the following properties: 

(1) Its half-life Ti/ 2 should be longer than its preparation 
and cooling time. 

(2) It should have a large fraction of the EC decay, im- 
plying large Z. 

(3) The life time r of the upper hyperfine state should 
be longer than T1/2, implying small Z. 

(4) The frequency f2 of the time modulation of the EC 
decay should be greater than the EC decay constant A. 
Table Q] shows H-like ions which may be considered as 
candidates for the experimental observation of time mod- 
ulation. The hyperfine splitting AE is calculated by Eq. 
([T]). The mean time r for the magnetic dipole transition 
between the upper and lower hyperfine states is estimated 
with the expression fl2j 

„2\ 3 

(15) 
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TABLE I: Hydrogen-like radioactive ions with allowed 
Gamow- Teller transition 1 + — ¥ + appropriate for observa- 
tion of the time modulation of the EC decay. The half-lives 
7i/2 of nuclear ground state is given for neutral atoms. 



H-like 


/' 


AE 


T 


21/2 


Qec 


EC//3+ 


ion 


fJ-N 


meV 






keV 


% 


18p 


+0.58" 


0.672 


5.55y 


109.8m 


1655 


3.3/96.7 


64 Cu 




8.95 


10.3h 


12. 7h 


1673 


43/18 


68 Ga 


0.0117 6 


0.596 


7.95y 


67.6m 


2921 


8.9/88.0 


78 Br 


0.13 b 


9.75 


15.9h 


6.46m 


3754 


8/92 


82 Rb 


+0.554 6 


48.0 


7.20m 


1.27m 


4378 


4/96 



"Calculated with Schmidt's factors and g". 

b The data from Ref. |13|| . No sign is given, if it is unknown. 



where A e is the electron Compton wavelength and F is 
the angular momentum of the ground state. The H-likc 
ion 68 Ga seems to be the best candidate. In Table |TT] 
we have compared the amplitudes of the two fields with 
different polarizations calculated from Eqs. (|10l) and (fT2")) 
for different photon number K and the fixed width ■ 
The trap field is Bq = 1.0 T. It is seen that the circularly 
polarized field is more suitable for experimental use. 

TABLE H: Strength of the linearly polarized (6; = b x for 
transition (i) with ujq /2n = 144 GHz) and circularly polarized 
(b c — b x = b z for transition (ii) with uio/2-n = 126 GHz) 
harmonic fields required to observe the A'-photon resonance 
with the half width r K /27r = 0.1 Hz in the ion 68 Ga 30+ . 



K 


1 


2 3 


4 


5 


bi,T 


1.54-10" 11 


2.38- 10" 3 




9.07-10 -2 


b c ,T 


1.00-1Q- 12 


2.62-10" 6 3.47-10" 4 


3.94-10" 


* 1.69-10 -2 



Storage ring ESR. — We try now to use multiple transi- 
tion mechanism for the interpretation of the GSI oscilla- 
tions. It should be emphasized, however, that this prob- 
lem is a more complicated one owing to the extremely 
complex configuration of the ring magnetic field designed 
to store ions. Moreover, the available experimental infor- 
mation is yet incomplete. 

The ESR lattice has periodic structure with the spatial 
period L = C/2, where C = 108.36 m is the circumfer- 
ence of the ring. Thus, an ion circulating in the ESR is 
exposed to the periodic magnetic fields of bending dipole 
and quadrupole magnets. We consider the two fields, 
which are in the direction of the x-axis tangent to the 
main orbit and the z-axis perpendicular to its plane. In 
the reference frame of a moving ion they can be expressed 
by the Fourier series 
00 

b x{t) = b xn sin (2im-fvt/L) , 

n=l 

00 

B z {t) = B z0 + Y / b zn cos(2nn 1 vt/L), (16) 
n=l 
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where t is the proper time. For the ions in consideration 
the Lorentz factor is 7 = 1.43. The frequency u> n of n-th 
harmonic is equal to nu>, where lo/2h = jv/L = 5.65 
MHz is the frequency of the first harmonic and the inter- 
val between nearby harmonics. The latter is much greater 
that the width of the multiphoton resonance at the fre- 
quency ui n . Therefore, the resonances with different n do 
not overlap, and may be considered as independent. 

At first we consider the transition, induced by the fields 
of n-th harmonic, between the hyperfine states F = 1/2 
and F = 3/2 with the energy AE ~ 1 eV. In principle, 
the resonance transitions with the Rabi frequency of 0.1 
Hz is possible for very high n and small K. However, 
this is a fallacious result, because the mean life time of 
the upper F — 3/2 state, r = 0, 03 s, is much shorter 
than the modulation period. Hence, the "sterile" state 
can not be the source of observed oscillations. 

Let us next consider the transitions between the Zee- 
man components of the lower hyperfine state. For the 
modulation to be possible, the weak interaction operator 
of the EC decay has to have the different matrix ele- 
ments for the states with different M. We assume that 
this is the case. Then the two quasi-stationary states 
exposed to the fields of n-th harmonic evolve according 
to where a+ and a_ are the amplitudes of the lower 
|l/2,+l/2) and upper 1 1/2, — 1/2) states with energies 
e\ = h(uj + — i\/2) and 63 = h(co- — i\/2), respectively. 

Suppose the ion 14 0p r 58+ ( Qr i42p m 60+^ c i rcu i atmg m 

the ESR moves along the main orbit. Then it is ex- 
posed to the periodic magnetic field of dipole magnets 
only. One period L contains three identical dipoles ar- 
ranged symmetrically around the center of the central 
dipole magnet. The dipole magnet of length d = 6.5 m 
involves the strong bending field B(0, 0, Bo), B = 1, 197 
T, and the fringe field h(b x ,b v ,b z ). The latter is cal- 
culated by the code OPERA [lj]. It is localized in the 
small region Ax = 29 cm on the two sides of each dipole. 
The component b z decreases monotonically from Bq to 
zero. The component b x has a pulse shaped form with 
the maximal value bo = 36 mT, and the component b y is 
too small to be taken into account. The Fourier coeffi- 
cients in (|16[) are easily found for these components and 

the matrix elements in (|2|) arc U+ = U- = pf n /n, V-\ = 

V-+ = —qg n /n, where 

g^Bbo , _ gnBjB d 



g^BlB Q 

p= , q = 



and 



(17) 



is the Zeeman splitting w = w_ — The periodic 

in n functions /„ ~ g n < 1 depend on the lattice pa- 
rameters. The splitting vj/2tt — 5.65 GHz is well above 
the frequency of the first harmonic. This means that 
the harmonics with high n are involved in the resonance 
transitions. Hence, the driving field is weak and we can 
use the approximate solution of ([7]) . In Table III we show 
the half width and detuning for the X-photon resonance 
calculated according to the formulas 



r K = 



gffn j, ( 2pf n 

irhn K \nhu>n 



for the harmonics, which are closest to the resonance. 



TABLE III: Half width and detuning for the harmonics close 
to the if -photon resonance for the ion 140 p r 58 + or 142 Pm 60+ 
moving along the main orbit. Here K is integer value of o/u„. 
All frequencies are in the laboratory frame. 



n 


K 


La- (Hz) 


A (MHz) 


n 


K 


V K (Hz) 


A (MHz) 


143 


7 


1.48 • 10" 10 


-6.92 


250 


4 


24.4 


-2.97 


166 


6 


2.88 ■ 1(T 5 


12.9 


251 


4 


0.304 


-18.8 


167 


6 


5.22 • 10" 4 


-10.9 


333 


3 


18.9 


0.989 


168 


6 


0.328 


-34.6 


500 


2 


854 


-2.79 


201 


5 


0.782 


-22.7 


999 


1 


2420 


0.989 



The large values of A in Table IHlI may be attributed to 
the inaccurate trajectory being used. In a storage ring, 
periodic motion of an ion in the axial direction (rotation) 
is coupled with its periodic motion in the transverse di- 
rection (betatron oscillations). The latter leads to an ion 
being exposed to the magnetic field of quadrupole lenses 
which change Tk and A. However, for the resonance 
condition A ~ Tk to be realized an ion has to move 
along a particular, "resonance" orbit. The decays of ions 
stored on non-resonance orbits generate the background 
with the conventional exponential law. 

The plausible explanation of the observed oscillation 
amplitude is related to the coupling of the F = 1/2 and 
F = 3/2 states. The Hamiltonian corresponding to Eqs. 
([7)1 offers the coupling due to virtual transitions in the 
states F=l/2 and 3/2 with different M and nonzero pho- 
ton numbers. The transitions change the energy (like 
the Lamb shift) and the wave functions of the states 
1 1/2, ± 1/2} . However, these changes are small due to 
the weak driving field of the n-th harmonic and cannot 
explain the observed amplitude. 

It is not inconceivable also that the coupling of these 
states is effected by the static or periodic potential W 
associated with the ESR fields. The solution of the 
Schrodiger equation for the two hyperfine states gives 
in the limit W^j^^ 2 / AE <C 1 (the matrix elements W 
is assumed to be real) the new eigenfunctions 

2(WW)l'K 



cos • 



II I) 

I 2' 2/ 
II 

! 2' I' 



, tan 9 - 



hmAE 



With this basis, we find, following the method of 
the quasienergy function 



*(t) = (6l|l> + &2|2)£ 



-iKuj n t\ 



f (lj— -\-u)+ — i\)t 



(20) 



, A = zu - Kuj n , 



(18) 



satisfying the initial condition ^(O) = |1). Here the co- 
efficients 

L A 

61 =cos (flt/2)-i-^- sin (fii/2) , b 2 = -i— sin (fit/2) 

represent Rabi oscillations (slow motion). The function 
(I2TJ1) describes the evolution of the system and allows to 
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calculate the amplitude of the transition of the mother 
ion from the state lsi/ 2 ,F — 1/2 into the daughter nu- 
cleus. The probability of the EC decay is proportional to 
the square of the amplitude averaged over the period of 
n-th harmonic (fast motion) 



P, 



EC 



\M\ 



1 — sin( 



2ri „ fn 



(21) 



where \M\ 2 = |M + i| 2 = |M_i| 2 are the squared am- 
plitudes of EC decay calculated with the unperturbed 
states. Equation (l2~Tj) describes damped Rabi oscillations 
with the EC decay constant A = \ec and the amplitude 
sin6> ~ W 2 /(hmAE). Comparing the latter with ob- 
served amplitude 0.2, we find W ~ 0, 001 eV. The elec- 
tric quadrupole field required to achieve this value is of 
the order of atomic fields and cannot exist in the ESR. 
The required magnetic field B(0, B y ,B z ) should be of the 
order of 20 T, which is greater by a factor 20 than the 
ESR fields. 

Conclusion. — The effect of atomic electron on the EC 
decay of H-like ions is very important. We have shown 
for the first time that the EC-decay rates of such ions can 
be modulated by using the Rabi resonance method with 
single or multiphoton transitions. For a long time the 
multiphoton resonance has been subject of intense exper- 
imental and theoretical studies in different fields except 
nuclear physics. We have proposed the experiment with 
intermediate mass ions stored in a Penning trap which 
can demonstrate this phenomenon. 



In an attempt to understand the GSI oscillations, we 
have proposed an alternative mechanism involving reso- 
nance multiphoton transitions, induced by the periodic 
magnetic field of the ESR, between the magnetic sub- 
states of the ground F = 1/2 state of the radioactive 
ions 140 p r 58 + and 142 Pm 60+ . The mechanism includes 
the coupling of F = 1/2 and F = 3/2 hyperfine states to 
explain the oscillation amplitude. The important point of 
this scenario is that the predicted EC decay constant does 
not violate experimentally and theoretically established 
ratio \^ Hke /^EC~ Hke «3/2 0. We find that the pe- 
riodic magnetic field of the ESR can generate the oscilla- 
tions with the frequency of about 0.1 Hz in the transitions 
involving several photons. However, the field required to 
explain the observed amplitude should be many times 
higher than that of the ESR. Moreover, to tune to the 
resonance a stored ion should move along specific, "reso- 
nance" orbits. This cast some doubts upon the involve- 
ment of the hyperfine structure in the observed oscilla- 
tions. Therefore, it is necessary to perform the single-ion 
experiment with the ions 123 > 5 Cs 54+ having the ground 
state lsi/ 2 F = or with the He-like ions without a hy- 
perfine structure. From these experiments definite con- 
clusions may be drawn regarding competing hypotheses 
for the GSI oscillations. 
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